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Abstract 

The  ability  to  compute  the  maximum  area  on  the 
earth’s  surface  (footprint)  reachable  by  an  autonomous 
air  vehicle  can  be  useful  in  planning  for  the  vehicle’s 
safe  operations.  The  information  can  be  important 
when  the  vehicle  experiences  subsystem  failures  caus¬ 
ing  it  to  be  unable  to  maintain  its  nominal  perfor¬ 
mance.  In  this  paper,  we  present  a  method  to  cal¬ 
culate  the  footprint  of  a  reusable  launch  vehicle  that 
experiences  a  failure  in  one  or  more  of  its  aero-control 
surfaces.  During  a  control  effector  failure,  the  max¬ 
imum  attainable  moments  of  the  vehicle  are  reduced 
which  may  decrease  the  range  of  conditions  that  the 
vehicle  can  maintain  a  trimmed  condition.  Addition¬ 
ally,  the  lift  and  drag  characteristics  of  the  vehicle  can 
change  when  control  effectors  are  moved  to  off-nominal 
positions  to  correct  for  moment  imbalance  caused  by 
failures  or  damage.  As  a  result,  the  footprint  of  the  ve¬ 
hicle  is  reduced.  A  technique  for  calculating  the  avail¬ 
able  effectiveness  of  the  aero-control  surfaces  is  used 
in  conjunction  with  a  footprint  generation  algorithm 
to  include  the  effects  of  rotational  trim  on  the  vehicle 
footprint. 


1  Introduction 

Future  space  operation  vehicles  aim  to  achieve  afford¬ 
able  and  reliable  access  to  space  by  minimizing  the 
weight  and  number  of  subsystems  and  fluids  used  in 
the  vehicle.  This  effort  hais  produced  a  number  of  ad¬ 
vances  in  propulsion,  structures,  materials,  and  system 
guidance  and  control  for  reusable  launch  vehicles.  Be¬ 
cause  of  the  high  cost  penalty  for  additional  weight,  it 
is  desirable  for  autonomous  hypersonic  vehicles  to  have 
limited  hardware  redundancy.  As  a  result,  the  system 
reliability  and  safety  associated  with  hardware  redun¬ 
dancy  may  be  greatly  reduced.  Furthermore,  control 
effectors  are  normally  sized  to  be  maximally  used  dur¬ 


ing  nominal  flight  conditions  so  that  actuator  hard¬ 
ware  limits  are  encountered  when  one  of  the  actuators 
fail.  Consequently,  a  single  control  actuator  failure  can 
severely  affect  the  vehicle’s  performance  and  safety.  To 
compensate  for  the  reduction  in  system  hardware  ro¬ 
bustness  to  control  effector  failures,  attention  is  be¬ 
ing  focused  on  system  software  by  designing  guidance 
and  control  laws  that  are  adaptive  to  these  failures  [1], 
[2],  [3].  Control  actuator  failures,  therefore,  along  with 
their  effects  on  mission  performance  and  safety  are  con¬ 
sidered  in  the  whole  process  of  trajectory  planning  and 
re-targeting.  An  important  element  of  the  vehicle’s  safe 
operation  in  contingent  flying  conditions  is  the  ability 
to  compute,  in  real  time,  the  largest  reachable  area  on 
the  earth’s  surface  given  its  current  conditions,  as  well 
as  constraints  on  skin  temperature,  structural  load, 
and  achievable  aerodynamics.  The  reachable  area  is 
referred  to  as  the  footprint  of  the  vehicle.  The  foot¬ 
print  information  is  useful  in  the  event  of  a  system 
failure  and  the  flight  path  of  the  vehicle  needs  to  be 
altered.  In  this  paper,  we  apply  a  technique  based  on 
linear  optimization  with  hard  constraints  on  the  vari¬ 
ables  to  compute  the  maximum  moments  achievable  by 
the  vehicle’s  current  control  effectors  to  compute  the 
footprint  of  a  hypersonic  vehicle  under  nominal  and 
contingent  operating  conditions  [5],  [6]. 


2  Problem  Formulation 

A  typical  trajectory  of  a  hypersonic  vehicle  generally 
consists  of  five  distinct  flight  segments:  ascent,  orbit, 
reentry,  terminal  area  energy  management,  and  ap¬ 
proach  and  landing  .  Each  flight  segment  is  charac¬ 
terized  by  unique  system  control  settings.  Early  in  the 
ascent  phase,  for  example,  the  vehicle  experiences  low 
dynamic  pressure  and  relies  mainly  on  propulsive  con¬ 
trol  and  thrust  vectoring  to  maintain  its  attitude  and 
trajectory.  Power-pack-out  or  engine  failure  is  an  im¬ 
portant  control  failure  to  consider  during  ascent.  On 
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the  other  hand,  during  the  reentry  phase,  the  propul¬ 
sion  system  is  normally  ineictive  and  not  part  of  the  ac¬ 
tive  control  effector  suite.  The  aero-control  surfaces  are 
the  control  effectors  that  are  used  to  direct  the  vehicle 
during  this  phase.  This  research  concentrates  on  com¬ 
puting  the  effects  of  control  failures  during  the  reentry 
phase.  The  type  of  failures  considered  include  symmet¬ 
ric  effector  failure  where  surfaces  are  locked  in  a  par¬ 
ticular  position.  In  particular,  a  method  for  estimating 
the  effect  of  symmetric  failures  on  force  and  moment 
equilibria  and  the  maximum  trimmable  lift  coefficient 
is  presented.  The  vehicle  model  used  in  this  study  is 
that  of  the  X-33  shown  in  Figure  1,  which  was  to  be 
a  technology  demonstrator  for  a  single  stage  to  orbit 
spacecraft. 


Figure  1:  X-33  Reusable  Launch  Vehicle 


2.1  Method  for  Determining  the  Range  of 
TVimmable  Angles  of  Attack 
We  will  consider  trim  conditions  where  the  angular  ve¬ 
locity  vector  is  zero.  We  will  therefore  neglect  the  ef¬ 
fects  of  rotational  damping,  thus,  the  rolling  moment 
L,  pitching  moment  M,  and  yawing  moment  N  of  the 
vehicle  about  its  body  axes  are  functions  of  the  angle 
of  attack  a,  side  slip  0,  altitude  h,  Mach  number  M, 
and  control  deflection  S: 


L  = 

L{a,p,h,M,S) 

(1) 

M  = 

M{a,p,h,M,S) 

(2) 

N  = 

N{a,p,  h,M,6) 

(3) 

These  moments  can  be  further  separated  into  those 
generated  by  the  base  vehicle  configuration,  i.e.,  wing 
body,  propulsion,  (Lo,  Mo,  No),  and  those  generated 
by  a  suite  of  m  control  effectors  (Lj,  Ms,  Ns)'. 

L  =  Lo{a,/3,h,M)  +  Lg{a,p,h,M,S)  (4) 
M  =  Moia,p,h,M)  +  Ms{a,p,h,M,S)  (5) 
N  =  No{a,P,h,M)  +  Ng{a,p,h,M,5)  (6) 

where 

Ls{a,P,h,  M,  6)  =  T,^iLsi{a,p,h,M,[0---5i---0]) 


where  Ls^  (a,  p,  h,  M,  6i)  denotes  the  rolling  moment  in¬ 
crement  induced  by  the  control  surface  deflected  at 
Si.  The  vehicle’s  attitude  is  maintained  by  the  eiero- 
control  effectors  S  so  that  L  =  0,  M  =  0,  N  =  0 
in  Equations  4,  5,  6.  The  calculation  of  the  foot¬ 
print  of  the  vehicle  is  then  the  largest  region  reach¬ 
able  by  the  vehicle  while  still  maintaining  its  desired 
angular  attitude  with  L  —  0,  M  =  0,  iV  =  0.  Un¬ 
der  an  aero-control  effector  failure,  it  is  assumed  that 
the  configuration-based  angular  moments  Lo,  Mo,  No 
remained  unchanged  while  the  aero-control  effective¬ 
ness  Ls,  Ms,  Ns  are  reduced.  For  the  X-33  vehicle 
with  eight  aero-control  surfaces,  it  is  necessary  to  cal¬ 
culate  the  maximum  attainable  moments  achieved  by 
the  functioning  aero-control  effectors.  The  calculated 
maximum  attainable  moments  are  then  used  to  predict 
the  maximum  lift  coefficient  that  is  attainable  while 
maintaining  rotational  equilibrium.  Since  the  maxi¬ 
mum  attainable  moments  vary  along  the  vehicle’s  trar 
jectory,  the  net  moment  generating  capability  of  the 
control  effectors  must  be  calculated  along  the  vehicle 
flight  paths  that  are  used  to  estimate  the  vehicle  foot¬ 
print.  Figure  2  gives  an  conceptual  example  of  how 
the  attainable  moment  set  might  change  under  a  con¬ 
trol  failure.  The  maximum  attainable  moments  for  a 


^ _ _  L 


Figure  2;  Attainable  Moment  Sets  for  a  Linear  Plant 

suite  of  linear  effectors  can  be  calculated  by  consider¬ 
ing  the  largest  moment  producing  capabilities  of  a  set 
of  control  effectors  [4].  Since  the  moment  generated  by 
each  control  surfeice  is  a  linear  function  of  the  surface 
deflection,  the  largest  moments  are  found  by  exam¬ 
ining  the  combined  generated  moments  at  the  upper 
and  lower  limits  of  their  deflection  ranges.  The  geo¬ 
metric  envelope  of  the  maximum  attainable  moments 
are  then  represented  as  facets  in  the  moment  space  as 
seen  in  Figure  2.  For  a  hypersonic  vehicle  travelling 
in  a  dynamic  environment,  the  moments  generated  by 
individual  control  effectors  are  often  nonlinear  func¬ 
tions  of  control  deflection.  The  method  presented  be¬ 
low  considers  only  symmetric  control  effector  failures. 
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This  condition  corresponds  to  cases  where  the  left  and 
right  control  surfaces  fail  in  an  identical  manner  so  that 
L  =  0  and  N  =  0.  Such  a  condition  is  not  as  restrictive 
as  it  may  seem  since  reconfigurable  control  laws  often 
balance  the  lateral  moments  caused  by  asymmetric  fail¬ 
ures  by  deflecting  the  opposing  control  surface  to  the 
same  (or  nearly  the  same)  deflection  as  the  failed  sur¬ 
face.  For  such  cases,  following  method  can  be  used  to 
calculate  the  attainable  moments  using  the  nonlinear 
aerodynamic  database  in  the  pitch  axis. 

First  note  that  in  Equation  5,  Ms  can  be  written  as 

Mg  =  \pv'^-cSCm{8) 

where  p  is  the  air  density,  v  is  the  velocity,  S  is  the  plane 
form  area,  c  is  the  mean  aerodynamic  chord,  Cm^  is  the 
total  control  pitching  coefficient.  The  control  pitching 
coefficient  can  be  further  expressed  as 

Cmq  =  C’mj,  H - f 

The  upper  bound  Cms  and  lower  bound  of  Cmc 

O  <)  w 

can  be  expressed  as 

Cm j  =  Cmi^  +  •  •  •  +  Cms^ 

Q-M^  =  +  •  •  •  + 

where 

Cm  Si  =  maxCm((5i) 

Cm  Si  =  minC„((5i) 

<>t 

Thus,  an  angle  of  attack  is  trimmable  if  Cmo  hes  be¬ 
tween  the  upper  bound  Cmx  and  lower  bound  C_Ms- 

a  o 

From  figure  3,  one  can  determine  that  the  range  of 
trimmable  angles  of  attack  of  the  X-33  for  a  flight  con¬ 
dition  of  Mach  12  with  the  failures  of  left  and  right 
flaps  fixed  at  10  degrees  .  The  base  pitching  moments 
CMo  is  bounded  by  the  upper  and  lower  bounds  of  Cg 
for  the  angles  of  attack  between  17  degrees  and  37  de¬ 
grees.  Thus,  the  set  of  trimmable  angles  of  attack  at 
is  defined  as  {q|17°  <  a  <  37°}.  Using  this  range  of 
trimmable  angles  of  attack,  we  can  then  find  the  con¬ 
dition  in  which  the  vehicle  is  trimmed  in  its  angular 
attitude,  balanced  with  lift  being  equal  to  its  weight, 
and  encounters  the  least  drag.  In  the  next  section,  we 
will  see  how  these  factors  affect  the  vehicle  reachable 
region  under  no  power. 

2.2  Footprint  Calculation 

Assuming  that  the  vehicle  angular  orientation  is  main¬ 
tained  by  its  inner  loop  attitude  controller,  the  motion 
of  the  un-powered  vehicle  over  a  non-rotating  earth  can 
then  be  modelled  as  a  point  mass; 

h  =  nsin(7)  (7) 


Figure  3:  X-33  Pitching  Moments  at  Mach  12 


= 

^  = 
V  — 

7  = 
^  = 


V  cos(7)  cos(V’) 

(i?o  +  h) 

V  cos(7)  sin('^) 

{Ro  -b  h) 

—D  /isin(7) 

m  {Ro  + 

L  cos((t)  p  cos(7) 
mv 

Z.sin((T) 


+ 


mv  cos(7) 


cos(7) 
v{Ro  +  h)^  '  (Ro  +  h) 
V  cos(7)  cos(V’)  tan((?!>) 
{Ro  +  h) 


(8) 

(9) 

(10) 

(11) 

(12) 


Where 


h  :  Altitude 
9  :  Longitude 
(j)  :  Latitude 
V  :  Velocity 
7  :  Flight  path  angle 

Ip  :  Heading  angle 
m  :  Mass  of  the  vehicle 
Ro  :  Earth’s  Radius 
p  :  Gravitational  Parameter 


The  total  energy  of  the  un-powered  vehicle  is  strictly 
decreasing  because  the  non-conservative  forces,  such  as 
the  aerodynamic  drag  forces  and  friction,  are  acting 
on  the  vehicle.  The  vehicle  footprint  then  consists  of 
points  on  the  earth’  surface  at  which  the  total  energy 
decreases  to  a  set  value.  When  its  energy  reaches  this 
value,  the  vehicle  then  enters  the  final  part  of  its  tra¬ 
jectory  called  terminal  area  energy  management  phase. 
Combining  its  velocity  v  and  altitude  h  in  an  energy- 
state  approximation  [7],  a  reduced  order  model  can  be 
obtained  to  simplify  the  vehicle  description.  The  spe¬ 
cific  energy  of  the  vehicle  can  be  expressed  as 


E=-v^  + 


(Ro  +  h)^ 


(13) 
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FVom  equations  7,  10,  and  13,  the  rate  of  decrease  in 
the  vehicle’s  energy  is 


dE  _  vD  2/xhusin(7) 
dt  m  {Ro  +  h)^ 


(14) 


Assuming  that  the  vertical  forces  acting  on  the  vehicles 
are  in  equilibrium  and  its  vertical  motion  frequency  dy¬ 
namics  are  much  faster  its  horizontal  motion  dynamics, 
the  flight  path  angle  can  be  taken  to  be  near  zero,  i.e., 
7  w  0.  Moreover,  since  the  acceleration  normal  to  the 
velocity  vector  is  small,  i.e.,  (u-y  »  0)  and  u  0,  the 
time-rate  of  change  of  the  flight  path  angle  can  then  be 
taken  to  be  essentially  zero,  i.e.,  7  »  0.  From  equation 
11,  we  have 


m  /  /X  ^ 

cos{(j)  \  (iJo  +  h)2  (iio  -\-h)} 


(15) 


With 


L  =  \pv‘^SaCL 


1 


D  =  ^fyv^SaCo,  and  Cd  -  +  kCl 


(16) 

(17) 


we  can  combine  equations  15  and  16  into  equation  17 
to  get 


D  =  + 


2m^kg^ 
pv^Sa  cos^  a 


(18) 


where  g  is  the  effective  gravity  constant: 


.  _  p _ 

Ro  "I"  k  Ro  “f*  h 

Given  the  initial  velocity  v{to)  =  Vo,  and  initial  alti¬ 
tude  h(to)  =  ho,  the  vehicle  speciflc  energy  E{t)  can 
be  calculated  according  equation  14; 


Eit)  =  Eo  +  dE=  ^  +  dE  (19) 

The  vehicle’s  velocity  v{t)  can  be  derived  from  the  cur¬ 
rent  specific  energy  E{t)\ 

(20) 


since  the  vehicle  flight  altitude  h  is  small  when  com¬ 
pared  to  the  earth’s  radius  Ro-  From  the  above  discus¬ 
sion,  the  unpowered  vehicle  under  energy-state  approx¬ 
imations  has  the  reduced-order  model  x  =  /(x,  u,  t)  of 
the  form: 


9  = 

V  cos(V') 

(21) 

(Ho  +  h)  cos{4>) 

= 

V  sin(^) 

Ro  "f*  hf 

(22) 

tP  = 

zg  V  cos(V»)  tan(<^) 

V  Ro  +  h 

(23) 

E  = 

—Dv 

(24) 

m 

where  x  =  [9,  <f>,  V’>  E]"^  is  the  vehicle  state  vec¬ 
tor.  The  control  inputs  u  =  \p  z]"^  are  the  air  den¬ 
sity,  p  =  p{h),  and  the  tangent  of  the  bank  angle, 
z  —  tan(<T).  The  optimization  objective  is  to  find  the 
control  vector  u  =  \p  z]  such  that,  at  time  t  =  tf, 
the  crossrange  position  specified  by  the  vehicle  lati¬ 
tude  is  maximized  for  a  given  downrange  value 
9f  =  eitf): 


The  initial  values  of  the  vehicle  states  are  taken  to  be 
0(to)  =  0,  <l>{to)  =  0,  ip{to)  =  0  and  E{to)  =  Eo- 

The  Hamiltonian  function  H  to  be  minimized  is 

H  =  G  +  X^f  (26) 

=  — +  Xg6  +  X,f,ij)  -b  XeE  (27) 

The  weighting  coefficients  A  are  the  Lagrangian  multi¬ 
pliers  to  be  determined.  The  necessary  conditions  for 
optimality  are 


dH 

dX 

(28) 

dx 

dt 

dH 

du 

=  0 

(29) 

Taking  variations  of  the  Hamiltonian  function  H,  en¬ 
forcing  the  condition  dH  =  0  at  the  stationary  point  of 
H,  and  the  Lagrangian  multipliers  do  not  vary  at  the 
boundary  produces 

[H5  t  -  X'^dftl  =  0 

The  final  heading  rpit/),  the  final  down-range  position 
9{tf)  are  unspecified.  The  initial  and  final  values  for 
the  Lagragian  multipliers  X{to)  are  chosen  to  be  zero. 
The  dynamics  of  X{t)  iov  to  <  t  <tf  can  be  found  from 
equations  28; 


dH 

Xe  = 

dH 

di>' 

dE 

dH 

\  _ 

dH 

d9  ’ 

= 

d<t> 

A,;,  =  — 


Combining  with  equation  24,  it  can  be  shown  that 


Xg{t)  — 

0 

(30) 

A^(f)  = 

1  —  cos{9f  —  9) 

(31) 

A^(t)  = 

—cos{<l>)  sin(0/  —  9) 

(32) 

XE{t)  = 

Hr  Xg9  -t-  X,p^  -f  X,j,‘4) 

E  ~  E 

(33) 

At  the  optimal  control  Uopt,  H{x,  Uopt,  A,  t)  =  0  and 

H{x,  Uopt,  X,  t)  <  H{x,  Uopt  +  Au,  A,  t)  (34) 
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that  is, 

i?,- (x, ,  A, t) -j-  X^Ei^XjUopti ^ 

Hr{x,  Uopt  +  A«,  A,  t)  +  \eE{x,  Uopt  +  Aw,  A,  t). 

(35) 

From  equations  33  and  34,  it  can  be  shown  that  Sub¬ 
stituting 

_  -^rC^i  '^opty  A,  t) 

E{x^  V>opty  t) 

into  equation  35,  we  have 

Hr  jXy  Uppf,  A,  t)  ^  HrjXyUopt  -b  Aw,  A,f)  ^  ^ 

E{x,UoptyXyt)  E{x,Uopt  +  Au,x,t)  ~ 

From  expression  36,  we  can  define  a  new  Hamiltonian 
function  that  relates  the  vehicle  states  x,  its  co-states 
A,  control  u  to  its  decreasing  rate  of  energy  E: 


_  Hr  _  -l-  Xg6  -f  X^^  +  Xy^ip 

~~E  ^ 


(37) 


Putting  values  for  the  states  6,  E  and  their  co¬ 
states  X,j,y  Afl,  Xyj,  into  equation  37,  we  have 


Hnew  — 


-vsin(tjj)  ,  » 
R«+h 


tisin(Vi) 
*  Ro+h 


lAJ-C  M  V  vcos(,tli)ta.n(it>) 

I  +'^n^\v{R„+h)^  R^+h  R„+h  J 


m  "I  JSa^  > 


(38) 

For  a  decreasing  total  energy  (£  <  0),  the  new  Hamil¬ 
tonian  Hnew  is  to  be  maximized  with  the  control  pa¬ 
rameters  being  the  air  density  p  and  the  bank  angle 
a.  For  a  hypersonic  vehicle  under  energy-state  approx¬ 
imations,  methods  suggested  by  Schultz  [5]  and  Vinh 
[6]  are  applied.  Since  the  control  parameter  p  does  not 
appear  in  the  numerator  of  Hnewy  Hnew  is  maximized 
when  the  air  density  p  is  minimized.  From  the  drag 
equation  17,  we  have 


dp 


—  0  =>  Popt  — 


V^SaCdo 


(39) 


The  optimum  bank  angle  cr  associated  with  z  =  tan(<7) 
is 


dH  _  Q  .  _  C03(^)  sin(9)(-^-^)^ 

cos(0)  sin(V')— sin(<^)  sin(ff)  cos(^) 

(40) 

with  0  =  Of  —  9,  R  =  Ro  +  h.  From  equations  39 
and  40  the  final  downrange  Of  =  9{tf)  of  the  vehicle 
is  needed  a  priori  to  generate  an  optimal  bank  angle 
<j)opt  and  altitude  popt  —  p{h)  to  achieve  the  largest  fi¬ 
nal  cross-range  </>/  =  (j>[tf).  To  generate  the  vehicle 
entire  footprint,  it  is  suggested  in  [6]  to  use  the  vehi¬ 
cle  heading  ipo  =  '’/’(to)  as  the  sweeping  parameters; 
180°  <  V’o  <  90°.  The  intermediate  downrange  6'f 


and  cross-range  4>'f  are  transformed  by  the  coordinate 
transformation  to  obtain  the  final  6j  and  </>/: 


tan(0/) 

sin(^/) 


tan(0})  cos(V’o)  -I- 


t&n{<t>'f)  s[n{ip'f) 


cos{6f) 

=  sm{(p'f)  cos{'ipg)  -  sm{6'f)  cos{4>'f)  sm{tp'f) 


Equations  39  and  40  give  the  optimal  cross-range  </'(t/) 
for  a  given  downrange  6/.  Iterations  on  the  initial  val¬ 
ues  of  Of  may  be  necessary  so  that  the  initial  guess 
matches  the  final  value  0{tf)  resulted  from  the  equa¬ 
tions  24. 


3  Approximate  Solution 

Online  footprint  computation  using  the  above  method 
may  be  slow  because  iterations  on  Of  are  necessary. 
An  approximate  solution  to  the  maximum  crossrange 
calculation  was  proposed  by  Vinh  [6]  for  the  optimal 
bank  angle  command  signal. 

J^oppro*  =  tan  ^  arctan(  ^  (41) 


where 


1 

10,/rcz 


The  advantage  of  this  approximate  solution  is  that  it 
does  not  require  an  initial  guess  of  0{tf)  and  subsequent 
iterations.  This  bank  angle  control  was  demonstrated 
in  [6]  with  good  results. 


4  Application  to  a  Reusable  Launch  Vehicle 

In  this  section,  we  apply  the  methods  for  the  foot¬ 
print  calculation  to  the  X-33  reusable  launch  vehicle. 
The  X-33  is  an  autonomous,  reusable  launch  vehicle 
that  has  two  linear  aerospike  engines  and  eight  aerody¬ 
namic  control  surfaces:  inner/outer  elevens,  rudders, 
flaps  as  shown  in  Figure  1.  At  the  beginning  of  the 
reentry  phase,  the  2500-slug  vehicle  attains  the  veloc¬ 
ity  of  10, 000  feet/second  at  an  altitude  of  180, 000  feet. 
With  this  initial  energy,  the  unpowered  vehicle’s  cross- 
range  is  calculated  until  its  final  energy  is  equivalent 
to  the  energy  at  the  final  speed  of  2500  feet/second 
and  50, 000  feet  of  altitude.  The  air  density  constraint 
which  is  related  to  thermal  constraint  on  the  vehicle 
is  p  <  7.397  X  10~®slugs/ft^.  Similarly,  structure- 
related  air  density  constraint  has  an  upper  limit  of 
3.981  X  10"®slugs/ft^. 

An  important  requirement  that  must  be  satisfied  in  cal¬ 
culating  the  vehicle  footprint  is  the  maintenance  of  lift 
to  effective-weight  equilibrium  (L=W)  while  banking 
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the  vehicle: 


L  = 


\p'SaCL  =  W  = 


mg 

cos((t) 


where  Sa  is  vehicle  planform  area  and  p  is  the  air  den¬ 
sity.  The  respective  normal  and  axial  force  coefficients 
Cat  and  Cj\  obtained  from  the  vehicle’s  aerodynamic 
table  are  transformed  into  the  corresponding  lift  and 
drag  coefficients  Cl  and  Co- 


Cl  =  cos(a)C/v  —  sin(a)Cyi 
Cd  =  sin(a)C'jv  -I-  cos(a)C^ 

A  root  solver  based  on  the  Secant  method  is  then  used 
to  minimize  the  residual  of  L-W^.  As  discussed  earlier, 
the  new  Hamiltonian  Hnew  in  Equation  37  is  optimized 
when  the  vehicle’s  drag  is  minimized.  Scanning  the  al¬ 
titude  that  is  20, 000  feet  above  and  below  the  vehicle’s 
current  altitude,  we  look  for  the  next  optimal  altitude 
command  that  minimizes  the  vehicle’s  drag  subject  to 
the  constraint  L  =  W.  Once  the  optimal  altitude  is 
found,  the  bank  angle  command  (t  is  calculated  using 
Vinh’s  control  law  in  Equation  41.  As  shown  in  Equa¬ 
tion  41,  the  bank  angle  command  is  a  function  of  the 
vehicle’s  current  parasitic  drag  Cd„  coefficient  and  in¬ 
duced  drag  parameter  k.  The  parasitic  drag  coefficient 
Cd„  in  Equation  17  for  the  X-33  vehicle  is  found  from 
the  aerodynamic  data  table  by  iterating  over  the  vehi¬ 
cle’s  angle  of  attack  a  at  a  given  Mach  number  to  find 
the  drag  coefficient  at  zero  lift: 

~  ^D<,lo£,(o,Moeh)=0 

A  third-order  polynomial  expression  is  then  used  to 
parameterize  Cdo  in  terms  of  the  vehicle’s  velocity  to 
give  the  following  Equation: 

Cd„  =  -2.8809  X  10-‘‘Mach^  -t-  7.4321  x  IQ-^Mach^ 
-6.0914  X  lO-^Mach  +  2.9697  x  lO”' 


The  parameter  k  in  Equation  17  can  be  found  by  look¬ 
ing  up  the  values  of  the  lift  and  drag  coefficients  from 
the  vehicle’s  aerodynamic  data  table  and  solving: 

Cp-Cp^ 

Cl 


5  Rotational  Equilibrium  Under  Failures 

In  addition  to  ensuring  that  the  lift  on  the  vehicle 
equals  its  weight,  rotational  equilibrium  must  also  be 
enforced  to  maintain  the  vehicle’s  attitude.  A  trim  rou¬ 
tine  is  used  to  find  the  aero-control  positions  that  are 
necessary  to  balance  the  base  pitching  moment  pro¬ 
duced  by  the  wing-body  portion  of  the  vehicle.  The 
lateral  directional  moments  resulting  from  the  wing- 
body  portion  of  the  vehicle  are  assumed  to  be  zero  ( 


i.e.  an  assumption  of  zero  steady-state  sideslip  and 
wing-body  symmetry).  In  the  trim  routine,  the  roll, 
pitch  and  yaw  control  effectiveness  of  each  aero-control 
surface  at  a  fixed  Mach  number  and  angle  of  attack  is 
found  from  the  aerodynamic  table  using  small  pertur¬ 
bations.  With  B  being  the  pitch  control  effectiveness 
matrix,  S  the  aero-control  deflections  and  Mo  the  base 
pitching  moment,  a  linear  programming  formulation  [3] 
is  used  to  find  S  such  that: 

min  J  =  \\B5  —  Mo\\i  (42) 

6 

subject  to 

where  8  and  5  are  vectors  that  represent  the  lower  and 
upper  limits  of  the  control  effectors.  Control  surface 
failures  such  as  locked  control  surfaces  are  acconuno- 
dated  by  setting  the  upper  and  lower  limits  for  the 
locked  effectors  equal  to  one  another.  By  solving  the 
minimization  problem  posed  in  Equation  42,  the  po¬ 
tentially  undesirable  moments  produced  by  locked  ef¬ 
fectors  are  automatically  taken  out  by  the  un-failed 
surfaces  whenever  it  is  physically  possible  to  do  so. 
One  additional  comment  regarding  Equation  42  is  in 
order.  The  linear  programming  problem  that  is  posed, 
assumes  that  the  moments  are  linearly  related  to  sur¬ 
face  deflections.  This  is  rarely  the  case  in  practice  and 
at  best,  at  a  particular  flight  condition,  the  moments 
are  nonlinear  functions  of  surfeice  deflection  and  in  or¬ 
der  to  find  the  deflections  that  produce  a  desired  mo¬ 
ment  one  must  find  8  that  solves  the  following  equation 
subject  to  8<  8  <  8 

Mo  =  f{8)  (43) 

In  general,  if  the  solution  to  Equation  42  which  we  will 
call  8ip  is  substituted  for  8  into  Equation  43  one  will 
find  that: 

Mo  ^  mp)  (44) 

An  iterative  approach  has  been  developed  here  that 
solves  multiple  linear  programming  problems.  The  ap¬ 
proach  finds  a  modified  value  of  base  pitching  moment, 
namely  Mo^^j  that  produces  a  control  deflection  vector 
^iPn,od  produces 

Mo  =  /(<5/p_J  (45) 

The  iterative  procedure  uses  the  following  update  rule 
to  compute  Mo„,a- 

+  uj{Mo  -  f{8ip^^,, ))  (46) 

where  k  indicates  the  A:th  iteration,  and  w  is  a  param¬ 
eter  that  affects  the  convergence  properties.  As  stated 
previously,  a  locked  control  surface  can  be  accommo¬ 
dated  directly  through  the  control  allocator  by  setting 
the  upper  and  lower  position  limits  equal  to  one  an¬ 
other.  This  ensures  that  the  remaining  control  surfaces 
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are  used  to  balance  the  possibly  undesirable  effects  of 
the  locked  surfaces  as  well  as  balance  the  wing-body 
moments.  In  most  cases,  one  or  more  locked  surfaces 
create  a  condition  that  requires  the  free  surfaces  to 
move  to  off-nominal  positions  to  counter  undesirable 
effects  of  the  locked  surfaces  to  maintain  rotational 
equilibrium.  This  ultimately  translates  into  force  per¬ 
turbations  that  increase  the  overall  drag  on  the  vehicle 
and  therefore  reduces  the  size  of  the  vehicle  footprint. 


6  Example;  Locked  Surfaces 

The  X-33  aerodynamic  model  is  used  in  this  example 
to  optain  representative  force  and  moment  data  for  an 
RLV.  The  case  under  consideration  uses  the  techniques 
described  above  to  compute  footprints  under  nominal 
conditions  and  in  the  case  of  right  inboard  and  out¬ 
board  elevons  both  locked  at  —30  degrees.  In  order  to 
maintain  rotational  equilibrium,  the  control  allocator 
automatically  deflects  the  left  inboard  and  outboard 
elevons  to  —30  degrees  to  balance  the  lateral  direc¬ 
tional  moments  and  deflects  the  body  flaps  to  balance 
the  pitching  moments  produced  by  the  other  control 
surfaces  and  the  wing-body.  This  trimming  operation 
is  performed  at  each  time  step  at  which  the  equations 
of  motion  and  control  commands  are  computed,  in  this 
case,  the  integration  time  is  1  second.  Figure  4  com¬ 
pares  the  footprint  of  the  nominal  and  failed  vehicle 
and  one  can  see  that  the  failed  vehicle  footprint  is  re¬ 
duced  in  size.  This  is  because  of  the  increase  in  drag 
due  to  the  off-nominal  control  surface  deflections  that 
are  required  to  maintain  rotational  equilibrium.  Figure 


Figure  4:  Footprint  Comparison  between  Nominal  and 
Failed  Conditions  with  Right  Inboard  and  Out¬ 
board  Elevons  Locked  at  —30° 

5  shows  the  commanded  altitude  that  yields  minimum 
drag  from  the  initial  energy  to  the  final  energy  to  yield 
the  largest  crossrange.  The  discrete  profile  of  com¬ 


manded  altitude  arises  from  the  increments  used  in  the 
altitude  sweep.  A  smoother  time  history  can  be  ob¬ 
tained  by  using  smaller  altitude  increments  or  using 
results  from  the  coarse  sweep  to  initialize  an  numeri¬ 
cal  optimizer  to  find  the  altitude  at  which  minimum 
drag  occurs.  Using  the  sweeping  method  to  initialize 
a  numerical  optimizer  reduces  the  likelihood  of  getting 
trapped  at  a  local  minimum.  Both  approaches  result 
in  increased  computational  time.  Figure  6  shows  the 
time  history  of  the  bank  angle  command  to  achieve  the 
maximum  crossrange.  The  entire  footprint  of  the  ve¬ 
hicle  is  then  generated  by  setting  the  initial  headiang 
V'(to)  from  90  degrees  to  0  degree  and  calculating  the 
vehicle  trajectory  given  the  initial  and  final  energies. 


Figure  5:  Altitude  Command  for  the  Max.  Crossrange 


Figure  6:  Bank  Angle  for  the  Max.  Crossrange 


7  Effect  of  Trimmable  Angle  of  Attack 

It  should  be  noted  that  control  surface  failures  can  also 
reduce  the  range  of  trimmable  angles  of  attack  on  some 
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vehicles.  This  can  cause  a  maximum  wing-body  lift  co¬ 
efficients  constraint  to  be  replaced  by  a  more  conser¬ 
vative  maximum  trimmable  lift  coefficient  constraint. 
Whether  such  constraints  become  active  or  not  is  de¬ 
pendent  upon  the  vehicle;  however,  loss  or  damage  to 
any  of  the  aero-control  surfaces  could  potentially  re¬ 
duce  the  maximum  trimmable  emgle  of  attack.  This  re¬ 
duction  can  decrease  the  meiximum  achievable  lift  that 
strongly  influences  the  cross-range  that  can  be  achieved 
by  banking  the  vehicle.  In  this  case  studied  here,  such 
a  constraint  was  never  active  because  the  X-33  requires 
less  control  power  to  trim  at  high  angles  of  attack  than 
at  low  angles  of  attack.  Additionally,  minimum  drag 
conditions  on  this  vehicle  generally  occur  at  high  alti¬ 
tude  and  at  high  angles  of  attack.  Nevertheless,  it  is 
interesting  to  look  at  the  reduction  in  the  size  of  the 
vehicle  footprint  that  would  result  from  a  hypothetical 
reduction  in  trimmable  angle  of  attack.  Such  a  case 
can  be  seen  in  Figure  7.  In  calculating  the  two  foot¬ 
prints  shown  in  Figure  7,  the  vehicle  initial  conditions 
remain  the  same  while  its  maximum  trimmable  angle 
of  attack  is  reduced  from  50  degrees  to  7  degrees. 
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Figure  7:  Footprint  Sensitivity  to  Varying  Constraints  on 
Max.  Angle  of  Attack 


8  Conclusion 


In  this  paper,  methods  to  calculate  the  largest  reach¬ 
able  area  are  presented.  The  optimal  method  is  based 
on  variational  calculus  and  requires  iterations  on  the 
initial  guess.  The  suboptimal  method  is  based  on  ap¬ 
proximations,  but  is  more  practical  for  implementation. 
A  numerical  example  is  used  to  show  how  the  vehicle’s 
largest  reachable  area  can  shrink  under  failure.  Such 
information  can  be  useful  in  the  process  of  selecting  a 
safe  and  available  landing  site  under  contingencies. 
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